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Abstract 

We explore different families of quasi-periodically Forced Logistic Maps for the existence 
of universality and self-similarity properties. In the bifurcation diagram of the Logistic Map 
it is well known that there exist parameter values s n where the 2 n -periodic orbit is su- 
perattracting. Moreover these parameter values lay between one period doubling and the 
next. Under quasi-periodic forcing, the superattracting periodic orbits give birth to two 
reducibility-loss bifurcations in the two dimensional parameter space of the Forced Logistic 
Map, both around the points s n . In the present work we study numerically the asymp- 
totic behavior of the slopes of these bifurcations with respect to n. This study evidences 
the existence of universality properties and self-similarity of the bifurcation diagram in the 
parameter space. 

Universality and self-similarity properties of uniparametric families of unimodal 
maps are a well known phenomenon. The paradigmatic example of this phenomenon 
is the Logistic Map l a (x) = ax(l — x). Given a typical one parametric family of 
unimodal maps {l a } a( ^i one observes numerically that there exists a sequence of 
parameter values {d n }„eN C I such that the attracting periodic orbit of the map 
undergoes a period doubling bifurcation. Between one period doubling and the 
next there exists a parameter value s n , for which the critical point of l Sn is a periodic 
orbit with period 2 n . One can also observe that 

lim ^~ C? "~ 1 = li m s n~ s n-l =6 = 4 6692Q (1) 

n^oo d n+ i - d n n^oo s n+ i - S n 

The convergence to this limit 5 (the so-called Feigenbaum constant) indicates a self- 
similarity on the parameter space of the family. On the other hand, the constant 
S is universal, in the sense that one obtains the same ratio S for any family of 
unimodal maps with a quadratic turning point having a cascade of period doubling 
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bifurcations. In this paper we explore if the same kind of phenomenon can be 
observed when the one dimensional maps is forced quasi-periodically. The answer is 
that universality and self-similarity do manifest, but they do it in different manners. 
Moreover, we show that they occur in a more "restrictive" class of maps, in the 
sense that the quasi-periodic forcing has to have a very particular form. We provide 
a theoretical explanation to this phenomenon in [20, 21, 22J (see also [19J ) . 

1 Introduction 

In the late seventies, Feigenbaum ([HE]) and Coullet and Treser ([21]) proposed at the same 
time the renormalization operator to explain the universal features observed in the cascade of 
period doubling bifurcations of the Logistic Map. This explanation was based on the existence 
of a hyperbolic fixed point of the operator with suitable properties. The first proof of the 
existence of this point and his hyperbolicity were obtained with numeric assistance [16} 15]. In 
[2"5] Sullivan generalized the operator and provided a theoretical proof of the hyperbolicity using 
complex dynamics. In [171 0] a nice summary of the one dimensional renormalization theory can 
be found, up to the date of their respective publication, as well as contributions to this theory. 

The kind of maps that we consider are maps in the cylinder where the dynamics on the periodic 
variable is given by a rigid rotation and the dynamics on the other variable is given by an 
unimodal one dimensional map plus a small perturbation which depends on both variables. 
These maps are known as quasi-periodically forced one dimensional maps and they have been 
extensively studied ( [101 EH EHl UH El El El E] ) with a focus on the existence of strange non- 
chaotic attractors. 

The paradigmatic example in our case of study is the Forced Logistic Map, which is the map on 
the cylinder TxK defined as 



where (a, e) are parameters and w is a fixed Diophantine number (typically it will be the golden 
mean). This family has interest not only for phenomena related with the existence of strange 
non-chaotic attractors ( [HI OS IB] ) but also as a toy-model for the truncation of period doubling 
bifurcation cascade ( |15U13j ). The term "truncation of the period doubling bifurcation cascade" 
refers to the fact that, when one fixes e and let a grow, the attracting set of the map undergoes 
only a finite number of period doubling bifurcations before exhibiting a chaotic behavior. This 
differs from the one dimensional case, where the number of period doubling bifurcations before 
chaos is infinite. 

In [13] we studied the truncation of the period doubling cascade for the map (0). We observed 
that the reducibility of the attracting set plays a crucial role. We computed bifurcation diagrams 
in terms of the dynamics of the attracting set, taking into account different properties, as the 
Lyapunov exponent and, in the case of having a periodic invariant curve, its period and its 
reducibility. One of these bifurcation diagrams is reproduced in Figure Q] (see Table [T] for the 
label of every color). 

Let d n be the parameter value where the attracting periodic orbit of the one dimensional map 
doubles from period 2 n to period 2 n+1 . Figure [1] reveals that from every parameter (a,e) = 
(d n , 0) of the map ((2|) it is born a period doubling bifurcation curve of the attracting set. Let s n be 
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2 2.2 2.4 2.6 28 3 3.2 3.4 3.6 

Figure 1: Diagram of the parameter space of the map ([2]) for u = ^ l ■ The axis correspond 
to the parameters a (horizontal) and e (vertical). For the correspondence of each color with the 
properties of the attractor see Table[TJ The black lines represent the period doubling bifurcations 
from period one to two (left) and two to four (right). 



Color 


Dynamics of the attractor 


Black 


Invariant curve with zero Lyapunov exponent 


Red 


Chaotic attractor 


Blue 


Non-chaotic non-reducible attractor 


Grey 


Non-chaotic reducible attractor 


White 


No attractor (divergence to — oo) 



Table 1: Color coding for the Figure [TJ 



the parameter value where the critical point of the (non- forced) one dimensional family is periodic 
with period 2 n . In the case of analytic maps in the cylinder, the reducibility loss of an invariant 
curve can be characterized as a bifurcation, for more details see definition 2.3 in |13j . Figure [1] 
also reveals that from every parameter value (a, e) = (s n ,0) two curves of reducibility loss are 
born. These curves correspond to a reducibility-loss and "reducibility-recovery" bifurcations of 
the 2 n -periodic invariant curve. In [20] we prove, under suitable hypothesis, that these curves 
exist. In Figure [1] we can observe that the period doubling bifurcation curve born around (d n , 0) 
is confined by one of the reducibility loss bifurcation curves born around (s n , 0) and another one 
born around (s n+ i,0). 

In this paper we approximate numerically slopes of the reducibility loss bifurcations curves 
described above and use them to explore for the existence of universality and self-similarity 
properties within the quasi-periodic one dimensional maps. A posteriori we know that the self- 
similarity properties exists between the bifurcation diagrams of the map ([2]) for different values of 
ui. For a better presentation of the concept, let us start the discussion looking for self-similarity 
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withing the bifurcation diagram of the map ([2]) with oj fixed. 



2 Description of the computations 



Consider the sequence of values s n for which the Logistic Map has a superstable periodic orbit as 
before. In Figure [1] we can observe how the two curves (namely S~ and S£) that are born at the 
parameter values (a,e) = (s n ,0) define a region of non-reducibility around this point. Assume 
that each of these curves correspond to reducibility-loss bifurcations and they can be written 
locally with a as a graph of e (in [2Qj we give concrete conditions under which this is true). In 
other words, there exist a neighborhood U n of s n , interval / = [0, p) and functions a : I — > R 
and (3 : I -)• R such that S~ n U = {(a n (e),e)| e G /} and 5+ n U = {(/3 n (s),s)\e G I}. 

Assume that the functions a n (e) and /3 n (e) are written at first order as 

«n(e) = s n + a n e + o(e), and /3 n (e) = s n + + o(e). 



Let s* be the limit of the parameter values s n . The self-similarity in the bifurcation diagram of 
the Logistic Map is manifested in the following way. One observes that s n — s* ~ 5(s n +i — s*). In 
terms of the parameter space this corresponds to the fact that the affine map L(a) = 5 (a— s*)+s* 
sends (approximately) the point s n to s n -\. We would like to find an analog self-similarity in 
the parameter space of the Forced Logistic Map ([2]). We can use the curves to detect this 
self-similarity. Therefore we look for an affine map of the kind 



e J \ 5i ) \ e ) \ 

such that it maps the curves S~ to S~_i (and respectively to S^i)- If we impose these 
conditions to the local parameterization (around Sj) of the curves considered above then we get 

a n _i + a^t + o(t) \ f S \ / a n + a' n s + o(s) - s* \ f s* 
t ) \ Sx J{ s ) + {0 

for any t and s. Then, replacing s = t/5\ in the first coordinate and equating terms in the order 
of t, we obtain 

a n = 5 a n -i + (1 - 8 )s* 

rJ — $°nr' 
"n — Si n+V 

Using these two equations we have that the value 5i can be estimated as 5\ ~ 5q f n . For 

a n-l 

So we recover the estimation So ~ a "~""~ 1 , which converges to the Feigenbaum constant d. 
Therefore we replace Sq by S in the estimation of 5%, obtaining S± ~ 5 J n . 



To obtain a numerical approximation of the values a' n we have computed points in the curve 
(a n (e),e) G S*~ for small values of s = 2~ k h n for different values of fc = 1,2, ...,M and h n a 
prescribed value depending on n (which have been decreased when we increased n). Then, we 
have used this set of points to estimate the value a' °"^ 2 h ^! and we did three steps 

Z fin 

of extrapolation to improve the accuracy of the results. We have also considered consecutive 
approximations to estimate the accuracy of the results. These computations have been done 
with quadruple precision (using the library [I]). 
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To compute each of the values (a n (e) , e) with good accuracy we have used the following pro- 
cedure. We have computed the 2 n periodic invariant curve of the system approximating it by 
its Fourier expansion. We use the same method described in Section 3.2.1 in |13| (originally 
from [3j) to continue the zero Lyapunov exponent but tuned to continue the reducibility loss 
bifurcation curve. To continue this curve we use the characterization of the reducibility loss 
bifurcation given by Definition 2.3 in [13j. 

3 Description of the results 

The values of a' n actually depend on to (the rotation number of the system). Hence, from now 
on, we write a' n = a' n (to). The estimated values of a' n (to) for the family (J2J when to = v ^~ 1 are 
shown in Table [21 We have included the ratios a' n (to)/a' n _ 1 (uj) in the third column of the table. 
The value e a in the fourth column corresponds to the estimated accuracy (in absolute terms) for 
the value of a' n . In Tables [3] and H] we included the same values for to = 2 v ^~ 1 and to = 4 V ^~ 1 . 

We have computed the approximate value of for the same values of to, but in all cases the 
value obtained has been equal (in the accuracy of the computations) to — a[. Because of that 
these results will be omitted of the discussion. 

We can observe in Tables [21 [3] and 2] that the ratios a' i (ui)/a' i _ 1 (to) do not converge to a constant. 
Therefore (a priory) it seems that there are no self-similarity of the bifurcation diagram. In fact, 
we will see in Section 13.21 that one has to do a more subtle analysis to uncover the self-similarity 
properties of the family. 

3.1 Universality of the ratio sequence 

A remarkable fact is that the ratios a' i (uj)/a' i _ 1 (to) on the third column of Table [2] are approx- 
imately the same values on the third column of Table [21 but shifted on the index value n by 
one position. Actually, the bigger is n the closer are the values. The same phenomenon can be 
observed in Tables [3] and [H 

Let us introduce some additional notation to follow with the analysis of this phenomenon. 
Consider F a quasi-periodic forced map as follows, 



where / is a C r map and to an irrational number. This map F can be identified with a pair 
(to, f) e Tn (R\Q) x C r (T xR,R). Consider a' n (to) the slope of the reducibility-loss bifurcation 
introduced before. Nothing ensures yet their existence for a general map F = (to,f), but 
whenever they exist they can be thought also depending, not only on to, but also on the function 
/ of the map F considered [22] ■ I* 1 other words we consider a' n (to) = a' n (to, /). 

Given two sequences {rj}j g z + and {si}i^i + of real numbers, we say that they are asymptoti- 
cally equivalent [26] if there exists < p < 1 and Kq such that 

\fi — ^ | < Kop 1 Vi € Z + . 

We denote this equivalence relation simply by Sj ~ rj. 



F : 



T x R -»• 



T x R 




(4) 
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Now we can resume the discussion on the phenomenon observed on Tables [3] and [H Our 
numerical computations suggest that 

a:_ 1 ( W) /)~«U(2 W ,/)' 1 > 

Consider two different values ujq and ui\ such that 2ujq = 2wi( mod 1) . If equation ([5]) is true, 
then a'j(wo, /)/a'j_ 1 (o;o, /) and a'^uii, f)/a' i _ 1 (uji, /) should be asymptotically equivalent. This 
can be checked numerically. In Table [5] we have recomputed the same values of Table [2] but this 
time for uj = ^ (and / the function associated to the map ([2]) as before). Again the results 
obtained show asymptotic equivalence. 

Given F a skew product map like (g]), let F 2 denote the map composed with itself. Concretely, 
the map F 2 is given as F 2 (6,x) = (9 + 2uj,f(9 + u, f(6,x))). To simplify the notation, let us 
denote by f 2 the map defined as f 2 (9,x) := f(9 + u,f(9,x)). Then we have that the map F 2 
corresponds to the pair (2a;, f 2 ). Moreover we have that the 2 k periodic curves are 2 k ~ l periodic 
curves of the map F 2 , therefore a'^ui, f) = a' i _ 1 (2uj, f 2 ). Hence the relation given by ([5]) can be 
rewritten as 

a{_ 1 (2 W ,/!»)~ a'^ujy [ ] 

Now we have that the relation given by ([5]) can be explained as a consequence of a much more 
general phenomenon. We believe that the sequence J^^tlf/) * s (asymptotically) universal, in 
the sense that it does not depend on the map /. This would imply ([6]) and consequently ([5]). 

To check the universality of the sequence, we consider the following map, 

9 = 9 + lo, 

x = ax(l — x) + e cos(27r6>). 

This map is like ([5} but with an additive forcing instead of a multiplicative one. In the literature, 
sometimes ([2]) is referred as the Driven Logistic Map and ([7|) is referred as the Forced Logistic 
Map. We do not do this distinction in this paper and we consider both as two different versions 
of the Forced Logistic Map. 

Note that both maps are in the class of quasi-periodically forced one dimensional unimodal maps, 
with a quasi-periodic forcing of the type h{x) cos(2ir9), where h is a function of one variable. 
This is certainly a very restrictive class of maps. In Section [3.31 we explore what happens when 
the quasi-periodic forcing is not of this form. 

We have computed the slopes a' n (u), the associated ratios a' n (uj) / a' n _ 1 (uj) and the estimation 
of the accuracy for the family ([7]) as we did before for the family ([2]). The results are shown in 
Table [6] for u = and Table [7] for u = 2^-. 

Now we can compare the sequences a'^uo) / a' i _ 1 (uJo) in Table [2] (respectively [3]) with the ones 
of Table [6] (respectively [7]). Again we can observe that both sequences have an equivalent 
asymptotic behavior for equal values of u. This agrees with the conjectured universal behavior. 

3.2 Self-similarity of the bifurcation diagram 

In order to find self-similarity properties of the parameter space we need to refine our analysis. 
Given a one dimensional map in the interval g : I — )■ J, its (doubling) renormalization is defined 
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n 


c4(w) 


<(^)/a'n-l(w) 


e a 





-2 


0000000000e+00 




- - - 


3 


Oe- 


15 


1 


-5 


8329149229e+00 


2 


.9164574614e+00 


1 


5e- 


14 


2 


-8 


4942599432e+00 


1 


.4562633015e+00 


8 


4e- 


13 


3 


-1 


6351279467e+01 


1 


.9249798777e+00 


7 


4e- 


15 


4 


-1 


1252460775e+01 


6 


.8817004793e-01 


3 


Oe- 


14 


5 


-1 


2243326651e+01 


1 


.0880577054e+00 


1 


6e- 


13 


6 


-1 


8079693906e+01 


1 


.4766978307e+00 


1 


6e- 


11 


7 


-3 


4735234067e+01 


1 


.9212291009e+00 


2 


Oe- 


12 


8 


-2 


9583312211e+01 


8 


.5168023205e-01 


2 


le- 


12 


9 


-4 


1569457725e+01 


1 


.4051657715e+00 


4 


2e- 


10 


10 


-7 


8965495522e+01 


1 


.8996036957e+00 


9 


le- 


11 


11 


-7 


4500733455e+01 


9 


.4345932945e-01 


8 


le- 


10 



Table 2: Approximate values of a' n {uj) for the map ([2]) for u = 



n 


a' n (u)) 


°4M/°4-iM 


e a 





-2 


0000000000e+00 






4 


le- 


16 


1 


-4 


7793787548e+00 


2 


.3896893774e+00 


9 


2e- 


14 


2 


-9 


9177338359e+00 


2 


.0751094117e+00 


8 


Oe- 


16 


3 


-6 


9333908531e+00 


6 


.9909023249e-01 


4 


3e- 


15 


4 


-7 


5678156188e+00 


1 


.0915028129e+00 


2 


4e- 


14 


5 


-1 


1183261803e+01 


1 


.4777397292e+00 


2 


3e- 


12 


6 


-2 


1488744556e+01 


1 


.9215095679e+00 


3 


Oe- 


13 


7 


-1 


8302110429e+01 


8 


.5170682641e-01 


3 


le- 


13 


8 


-2 


5717669657e+01 


1 


.4051750893e+00 


6 


le- 


11 


9 


-4 


8853450105e+01 


1 


.8996064090e+00 


2 


6e- 


11 


10 


-4 


6091257360e+01 


9 


.4345961772e-01 


1 


2e- 


10 


11 


-7 


1498516059e+01 


1 


.5512381339e+00 


4 


le- 


09 



Table 3: Approximate values of ot' n (uj) for the map ([2]) for uj = 2-^ 



n 


a' n (u)) 




e a 





-2 


0000000000e+00 






5 


3e- 


15 


1 


-6 


1135714539e+00 


3 


.0567857269e+00 


4 


0e- 


17 


2 


-4 


6432689399e+00 


7 


.5950186809e-01 


8 


8e- 


16 


3 


-5 


1662366620e+00 


1 


. 1126292121e+00 


5 


4e- 


15 


4 


-7 


6637702641e+00 


1 


.4834338350e+00 


5 


2e- 


13 


5 


-1 


4738755184e+01 


1 


.9231728870e+00 


6 


7e- 


14 


6 


-1 


2555429245e+01 


8 


.5186497017e-01 


6 


9e- 


14 


7 


-1 


7643286059e+01 


1 


.4052316105e+00 


1 


3e- 


11 


8 


-3 


3515585777e+01 


1 


.8996226477e+00 


5 


8e- 


12 


9 


-3 


1620659727e+01 


9 


.4346134772e-01 


2 


6e- 


11 


10 


-4 


9051192417e+01 


1 


.5512387420e+00 


9 


Oe- 


10 


11 


-9 


2911119039e+01 


1 


.8941663691e+00 


5 


6e- 


12 



Table 4: Approximate values of a' n {uj) for the map ([2]) for u = 
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n 


c4(w) 


«n(w)/a' n _i(w) 


e a 





-2 


0000000000e+00 




- - - 


5 


3e- 


17 


1 


-3 


7459592187e+00 


1 


.8729796094e+00 


3 


6e- 


15 


2 


-6 


1798012892e+00 


1 


.6497246575e+00 


2 


4e- 


13 


3 


-1 


2170152952e+01 


1 


.9693437350e+00 


2 


le- 


15 


4 


-8 


4095313813e+00 


6 


.9099635922e-01 


9 


5e- 


15 


5 


-9 


1576001570e+00 


1 


.0889548706e+00 


5 


2e- 


14 


6 


-1 


3525371011e+01 


1 


.4769558377e+00 


5 


Oe- 


12 


7 


-2 


5986306241e+01 


1 


.9213008072e+00 


6 


5e- 


13 


8 


-2 


2132200091e+01 


8 


.5168703416e-01 


6 


6e- 


13 


9 


-3 


1099463199e+01 


1 


.4051681745e+00 


1 


3e- 


10 


10 


-5 


9076676864e+01 


1 


.8996043914e+00 


5 


6e- 


11 


11 


-5 


5736446311e+01 


9 


.4345940345e-01 


2 


5e- 


10 



Table 5: Approximate values of a' n (uo) for the map ([2]) for to = 



n 


a' n (u)) 


a n(^)/«n-l( w ) 


e a 





-4 


0000000000e+00 






2 


7e- 


14 


1 


-8 


1607837043e+00 


2 


.0401959261e+00 


5 


le- 


14 


2 


-1 


1166652707e+01 


1 


.3683309241e+00 


2 


4e- 


12 


3 


-2 


1221554117e+01 


1 


.9004400578e+00 


2 


le- 


14 


4 


-1 


4564213015e+01 


6 


.8629342294e-01 


8 


6e- 


14 


5 


-1 


5837452605e+01 


1 


.0874224778e+00 


4 


6e- 


13 


6 


-2 


3384207858e+01 


1 


.4765132021e+00 


4 


5e- 


11 


7 


-4 


4925217655e+01 


1 


.9211776567e+00 


5 


8e- 


12 


8 


-3 


8261700375e+01 


8 


.5167534788e-01 


5 


9e- 


12 


9 


-5 


3763965691e+01 


1 


.4051640456e+00 


1 


le- 


09 


10 


-1 


0213020106e+02 


1 


.8996031960e+00 


1 


le- 


10 


11 


-9 


6355685578e+01 


9 


.4345927630e-01 


2 


3e- 


09 



Table 6: Approximate values of a' n for the map ([7]) for to = 



n 


a' n {(o) 




e a 





-4 


0000000000e+00 






3 


8e- 


15 


1 


-6 


2417012728e+00 


1 


.5604253182e+00 


2 


4e- 


13 


2 


-1 


2036825830e+01 


1 


.9284527253e+00 


2 


le- 


15 


3 


-8 


2891818940e+00 


6 


.8865180997e-01 


9 


Oe- 


15 


4 


-9 


0187641307e+00 


1 


.0880161934e+00 


4 


9e- 


14 


5 


-1 


3318271659e+01 


1 


.4767291245e+00 


4 


7e- 


12 


6 


-2 


5587438370e+01 


1 


.9212281462e+00 


6 


le- 


13 


7 


-2 


1792312360e+01 


8 


.5168011135e-01 


6 


2e- 


13 


8 


-3 


0621808757e+01 


1 


.4051656497e+00 


1 


2e- 


10 


9 


-5 


8169300479e+01 


1 


.8996036760e+00 


5 


2e- 


11 


10 


-5 


4880369083e+01 


9 


.4345932701e-01 


2 


3e- 


10 


11 


-8 


5132515708e+01 


1 


.5512380316e+00 


8 


2e- 


09 



Table 7: Approximate values of a' n for the map ([7]) for to = 2 
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n 












— <Jl,n-l 


1 


1 


36175279e+01 


- 


- - 


1 


11579415e+01 


- 


- - 


2 


8 


29844510e+00 


-5 


3e+00 


7 


57460662e+00 


-3 


6e+00 


3 


7 


69807112e+00 


-6 


0e-01 


6 


97211386e+00 


-6 


0e-01 


4 


7 


57782290e+00 


-1 


2e-01 


6 


83972864e+00 


-1 


3e-01 


5 


7 


55390503e+00 


-2 


4e-02 


6 


81347460e+00 


-2 


6e-02 


6 


7 


54857906e+00 


-5 


3e-03 


6 


80758174e+00 


-5 


9e-03 


7 


7 


54747726e+00 


-1 


le-03 


6 


80631795e+00 


-1 


3e-03 


8 


7 


54724159e+00 


-2 


4e-04 


6 


80604419e+00 


-2 


7e-04 


9 


7 


54719154e+00 


-5 


Oe-05 


6 


80598601e+00 


-5 


8e-05 


10 


7 


54718076e+00 


-1 


le-05 


6 


80597353e+00 


-1 


2e-05 


11 


7 


54717846e+00 


-2 


3e-06 


6 


80597086e+00 


-2 


7e-06 



Table 8: Estimations of the value 5i, n = 5\ in {bj) = 8 a " n ^2°J ) ^ or u ° = 2 1 on ^ ne ano - 
wi = 2^^- on the right. 



n 


#l,n(wo) 


8l,n 


- <5l,n-l 


1 


9.52608610e+00 






2 


8.35338804e+00 


-1 


2e+00 


3 


8.23204689e+00 


-1 


2e-01 


4 


8.20385506e+00 


-2 


8e-02 


5 


8.19937833e+00 


-4 


5e-03 


6 


8.19817945e+00 


-1 


2e-03 


7 


8.19796400e+00 


-2 


2e-04 


8 


8.19791815e+00 


-4 


6e-05 


9 


8.19790879e+00 


-9 


4e-06 


10 


8.19790672e+00 


-2 


le-06 


11 


8.19790628e+00 


-4 


4e-07 



Table 9: Estimations of the value <5i jTl = Si >n (u) = d J 1 '"^!^ for w o = for the family © 

with w = 2 1 ■ 



as TZ(f) = A^ 1 o go go A with A an affine transformation. Given a q.p. forced map like one 
might try to define a renormalization for these kind of maps. The most simple choice would be 
to define its (doubling) renormalization also as T(F) = A^ 1 o F o F o A, with A a suitable affine 
map. If the map F has rotation number uj, its (doubling) renormalization will have rotation 
number 2oj. This is due to the composition of F with itself when we define its renormalization. 
This argument become more clear in the rigorous definition of the renormalization operator for 
quasi-periodically maps done in [20] , 

Consider the reducibility loss curves S^; as in the parameter space of the map ([2]) (see Section [2]). 
We have that these curves depend on the rotation number (i.e. S 1 ^ = S^(u)). Now we can look 
for self-similarity of the parameter space as in Section [2] but taking into account the doubling of 
the rotation number. In other words we can look for affine relationship between the bifurcation 
diagram around s n for ojq and the bifurcation diagram around s n _i for 2wo- Consider an affine 
map like ([3]), but such that it maps the curves S~ (u) to S^_ 1 (2uj) (and respectively S^~(uj) to 
S' i t 1 (2a;)). Redoing the same computation of Section [2] we have that if there exists an affine 

relation between the parameters spaces it should be given by 5\ 8 tt '^J U . 
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0.03 



0.025 



0.02 



0.015 



0.01 



0.005 




3.25 3.3 3.35 3.4 3.45 3.5 3.55 




IA 2,6 



3,4 3.6 



0,0014 



0,0012 



0.001 



0,0008 



0,0006 



0,0004 



0,0002 




0,012 



0.01 



O.OOi: 



0,006 



0,004 



0,002 




3.5 3.51 3.52 3.53 3,54 3,55 3,56 3.57 



3.25 3.3 3.35 3.4 3.45 3.5 3.E 



Figure 2: Diagram of the parameter space of the map ([2]) for to = ^ on the hand side 
and for ui = — 1 on the right. See Table Q] for the coding of each color and the text for more 
details. 



In Table [8] we have estimations of these values for the family ([2]) for oj = ^ 1 (left) and 

u = \/5 — 1 (right). In Table [9] we have the same estimation for the family ([7]) and oj = v/ ^~ 1 . In 
all cases the sequence converges to a concrete value. This convergence means that there exists an 
affine relation between the reducibility loss bifurcation curves around the superstable periodic 
orbits of the uncoupled map. Note that the limit constants obtained in each of these three cases 
are different one to each other. This indicates that the renormalization factor depends both on 
the value of u) taken and the family of maps considered. Therefore the renormalization factor 
5\ is not universal. 

These evidences of self-similarity are only valid for infinitely small values of the coupling pa- 
rameter e. If the self-similarity between families extends to larger values of e, the bifurcation 
diagram of the map ([2]) in a box I\ x I2 of the parameters space for a prescribed U) = ujq should 
be approximately the same as the diagram in the box L(I\ x I2) for uj = 2ujq, where L is the 
affine map ([3]). 

In Figure [2] we have a bifurcation diagram of the map ^ analog to the one displayed in Figure 
[TJ The boxes have been selected such that the one in the left is the image of the one in the right 
through the affine map L given by ([3]). The value of <5o has been taken equal to the Feigenbaum 
constant (~ 4,66920) and 5± ~ 7,54718 the experimental value obtained in Table [HJ The 
results indicate that self-similarity properties extend to the whole reducibility region around 
each period doubling bifurcation. 
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3.3 Non-universality of the ratio sequence under more general conditions 

A natural question to ask after the numerical evidences of universality and renormalization 
reported in the previous sections is "how general are these phenomena?". In this section we 
present an example which demonstrates that the universality and the self similarity properties 
depend on the Fourier expansion of the quasi-periodic forcing. We must say that we designed 
this example after developing some of the theory presented in [201 12T1 [22] and the cited theory 
provides a theoretical explanation to both behaviors. 

Consider the following family, 

e = e + u, 

x = ax(l — x) + e(cos(2ir9) + E cos(47T0)). 

We consider E a fixed value and a and e as true parameters, obtaining a two parametric family. 
Remark that for E = we recover the family ([7J introduced in Section 13.11 

We have done the computation of the values a' n (cu) for the family (JSj) ? for E = 1CP 1 , 1CP 2 , 1CP 3 
and u> = 2^jfi. The results are shown in Tables [10] to [15] To compute the values in the 

tables we have used the same procedure that we have used for the families ([2]) and ([7]) before. 
In these tables we have also included the estimated values of the ratios a' n (u}) / a' n _ 1 (u;) and 
the estimated accuracies. In Table [TU1 we include the ratios 8i >n {oS) = a' n (uj) / a' n _ 1 (2u)) and the 
differences 5\^ n {uj) — 8i in -\{oS). 

In the third column of Tables [TU1 to [TBI we can observe that the sequence of ratios a' n (uj) / a' n _ 1 (oj) 
ceases to be universal. In other words, the sequence is not asymptotically equivalent to the 
sequences obtained for the maps ([2]) and ([7]) (displayed in Tables [2] to [7J . We can also observe 
in Table [16] that the different sequences Sa' n (u})/a' n _ 1 (2u}) cease to converge. Recall that the 
limit of this sequence gives us the scale factor between the bifurcations diagram of the map and 
itself for a doubled period. In other words, the self-similarity properties of the maps disappear. 

Analyzing the results with more detail we can observe in Tables [10] to [16] that (when the param- 
eter E is small) the map is not self-similar, but it behaves "close to self-similar" in the following 
sense. The values a' n (uj) of the family ([8]) (see Tables [TU1 to I15p differ form the same values of the 
family (J7J) (see Tables [2] and [3]) an order of magnitude similar to the magnitude of E. The same 
happens with the sequences a' n (uj) / a' n _ 1 (u) and a' n (uj)/a' n _ 1 (2u). Actually, this indicates that 
if two maps are "close", they still being "close" after several renormalizations (or after several 
magnifications of the parameter space). 

4 Summary, conclusions and further development 

In this paper we have done a numerical study of the asymptotic behavior of the slopes a' n (u>, f) of 
the reducibility loss bifurcations of quasi-periodic perturbations of the Logistic Map. Concretely 
we have considered families of maps in the cylinder TxR which can be written as, 

e = 6 + lu, 

x = ax(l — x) + eg(6, x), 
with cj an irrational number. Our numerical discoveries can be summarized as follows. 
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n 


a' n (uj) 


a n( w )/«n-l( w ) 


e a 





-4 


4000000000e+00 




- - - 


2 


le- 


14 


1 


-8 


5708073961e+00 


1 


.9479107718e+00 


3 


2e- 


14 


2 


-1 


2367363641e+01 


1 


.4429636637e+00 


2 


5e- 


12 


3 


-2 


2002414361e+01 


1 


.7790707058e+00 


2 


Oe- 


14 


4 


-1 


5466051366e+01 


7 


.0292519321e-01 


1 


2e- 


13 


5 


-1 


7124001858e+01 


1 


. 1071993396e+00 


7 


5e- 


13 


6 


-2 


5233583736e+01 


1 


.4735798293e+00 


5 


Oe- 


11 


7 


-4 


5526415150e+01 


1 


.8041993411e+00 


4 


4e- 


12 


8 


-4 


0793050977e+01 


8 


.9603037802e-01 


3 


7e- 


12 


9 


-5 


9579098646e+01 


1 


.4605207804e+00 


1 


le- 


09 


10 


-1 


0695246126e+02 


1 


.7951339260e+00 


9 


Oe- 


11 


11 


-1 


0464907069e+02 


9 


.7846341692e-01 


3 


le- 


09 



Table 10: Approximate values of a' n (uj) of the family ([8]) for oj = 1 and E = 10 1 . 



n 


a' n (u) 




e a 





-4 


4000000000e+00 






4 


2e- 


15 


1 


-6 


8849393233e+00 


1 


.5647589371e+00 


2 


6e- 


13 


2 


-1 


2339134533e+01 


1 


.7921921972e+00 


1 


9e- 


15 


3 


-8 


8252833730e+00 


7 


. 1522709709e-01 


1 


Oe- 


14 


4 


-9 


9507316299e+00 


1 


. 1275254526e+00 


2 


6e- 


14 


5 


-1 


5240261120e+01 


1 


.5315719172e+00 


4 


9e- 


12 


6 


-2 


6987715314e+01 


1 


.7708171207e+00 


6 


9e- 


13 


7 


-2 


4101506278e+01 


8 


.9305471016e-01 


1 


Oe- 


12 


8 


-3 


4274713216e+01 


1 


.4220983876e+00 


1 


4e- 


10 


9 


-5 


9739968275e+01 


1 


.7429750002e+00 


4 


Oe- 


11 


10 


-6 


0772668123e+01 


1 


.0172865818e+00 


3 


5e- 


10 


11 


-9 


7063373626e+01 


1 


.5971550472e+00 


8 


3e- 


09 



Table 11: Approximate values of a' n (oj) of the family (JSj) for oj = 1 , 2 and E = 10 1 . 



n 


a' n {uj) 


«n(w)/<-l(w) 


e a 





-4 


0400000000e+00 






2 


6e- 


14 


1 


-8 


1970065912e+00 


2 


.0289620275e+00 


4 


9e- 


14 


2 


-1 


1286613379e+01 


1 


.3769189098e+00 


2 


4e- 


12 


3 


-2 


1298993542e+01 


1 


.8871022535e+00 


2 


le- 


14 


4 


-1 


4654396022e+01 


6 


.8803232383e-01 


8 


9e- 


14 


5 


-1 


5964802815e+01 


1 


.0894207302e+00 


5 


Oe- 


13 


6 


-2 


3545880397e+01 


1 


.4748619617e+00 


4 


5e- 


11 


7 


-4 


4967877741e+01 


1 


.9097981041e+00 


5 


6e- 


12 


8 


-3 


8501241787e+01 


8 


.5619432630e-01 


4 


8e- 


12 


9 


-5 


434541 1904e+01 


1 


.4115236128e+00 


1 


le- 


09 


10 


-1 


0260055779e+02 


1 


.8879341272e+00 


1 


Oe- 


10 


11 


-9 


7178352178e+01 


9 


.4715227943e-01 


2 


3e- 


09 



Table 12: Approximate values of a' n (uj) of the family ([5D for uj = ^\ 1 and E = 10 2 . 



12 



n 


a' n (uj) 




e a 





-4 


0400000000e+00 




- - - 


3 


8e- 


15 


1 


-6 


2967793362e+00 


1 


.5586087466e+00 


2 


4e- 


13 


2 


-1 


2062745707e+01 


1 


.9157008786e+00 


2 


le- 


15 


3 


-8 


3390918858e+00 


6 


.9130959805e-01 


9 


Oe- 


15 


4 


-9 


1092083960e+00 


1 


.0923501648e+00 


4 


6e- 


14 


5 


-1 


3509885744e+01 


1 


.4831020608e+00 


4 


7e- 


12 


6 


-2 


5722289496e+01 


1 


.9039605503e+00 


6 


le- 


13 


7 


-2 


2023032848e+01 


8 


.5618478288e-01 


6 


8e- 


13 


8 


-3 


0953588852e+01 


1 


.4055098163e+00 


1 


2e- 


10 


9 


-5 


8281365651e+01 


1 


.8828629510e+00 


5 


le- 


11 


10 


-5 


5448581611e+01 


9 


.5139468666e-01 


2 


5e- 


10 


11 


-8 


6313143445e+01 


1 


.5566339289e+00 


8 


2e- 


09 



Table 13: Approximate values of a' n {ui) of the family ([8]) for uj = 2 ^ 1 and E 



n 


a' n {u) 




e a 





-4 


0040000000e+00 






2 


7e- 


14 


1 


-8 


1643491058e+00 


2 


.0390482282e+00 


5 


Oe- 


14 


2 


-1 


1178647202e+01 


1 


.3692024995e+00 


2 


4e- 


12 


3 


-2 


1229290655e+01 


1 


.8990930004e+00 


2 


le- 


14 


4 


-1 


4573231305e+01 


6 


.8646812284e-01 


8 


6e- 


14 


5 


-1 


5850170411e+01 


1 


.0876222355e+00 


4 


7e- 


13 


6 


-2 


3400073191e+01 


1 


.4763294391e+00 


4 


5e- 


11 


7 


-4 


4929300750e+01 


1 


.9200495820e+00 


5 


7e- 


12 


8 


-3 


8285483755e+01 


8 


.5212730036e-01 


5 


8e- 


12 


9 


-5 


3822109356e+01 


1 


.4058098286e+00 


1 


le- 


09 


10 


-1 


0217709608e+02 


1 


.8984223641e+00 


1 


Oe- 


10 


11 


-9 


6437862923e+01 


9 


.4383053171e-01 


2 


2e- 


09 



Table 14: Approximate values of a' n (oj) of the family (JSj) for uj = ^ 1 and E 



n 


a' n {uj) 


«n(w)/«'«-l( u ) 


e a 





-4 


0040000000e+00 






3 


7e- 


15 


1 


-6 


2470810635e+00 


1 


.5602100558e+00 


2 


4e- 


13 


2 


-1 


2039370639e+01 


1 


.9271993618e+00 


2 


le- 


15 


3 


-8 


2941266770e+00 


6 


.8891696467e-01 


9 


Oe- 


15 


4 


-9 


0277699390e+00 


1 


.0884533466e+00 


4 


8e- 


14 


5 


-1 


3337423899e+01 


1 


.4773774686e+00 


4 


7e- 


12 


6 


-2 


5600860673e+01 


1 


.9194756698e+00 


6 


le- 


13 


7 


-2 


1815381589e+01 


8 


.5213469453e-01 


6 


3e- 


13 


8 


-3 


0654499563e+01 


1 


.4051782426e+00 


1 


2e- 


10 


9 


-5 


8180013759e+01 


1 


.8979273708e+00 


5 


2e- 


11 


10 


-5 


4936892373e+01 


9 


.4425712239e-01 


2 


4e- 


10 


11 


-8 


5250385421e+01 


1 


.5517875464e+00 


8 


2e- 


09 



Table 15: Approximate values of a' n (oj) of the family ([8]) for uj = 1 \ and E 



13 



n 


5 


l,n(^o) 




— <5l,n-l 


1 


9. 


095188 




- - 


2 


8. 


387251 


-7 


.le-01 


3 


8. 


325844 


-6 


. le-02 


4 


8. 


182639 


-1 


.4e-01 


5 


8. 


035129 


-1 


.5e-01 


6 


7. 


730884 


-3 


.Oe-01 


7 


7. 


876621 


1. 


5e-01 


8 


7. 


902866 


2. 


6e-02 


9 


8. 


116387 


2. 


le-01 


10 


8. 


359271 


2. 


4e-01 


11 


8. 


040253 


-3 


.2e-01 



n 




<5l,n 


— <5l,n-l 


1 


9.473633 




- - 


2 


8.369274 


-1 


. le+00 


3 


8.244333 


-1 


.2e-01 


4 


8.205249 


-3 


.9e-02 


5 


8.183245 


-2 


.2e-02 


6 


8.137779 


-4 


.5e-02 


7 


8.162729 


2. 


5e-02 


8 


8.162820 


9. 


le-05 


9 


8.197747 


3. 


5e-02 


10 


8.219826 


2. 


2e-02 


11 


8.183173 


-3 


.7e-02 



n 


5 


l,n(<^o) 


<5l,n 


— <5l,n-l 


1 


9. 


520727 




- - 


2 


8. 


355159 


-1 


.2e+00 


3 


8. 


233307 


-1 


.2e-01 


4 


8. 


204041 


-2 


.9e-02 


5 


8. 


197777 


-6 


.3e-03 


6 


8. 


191961 


-5 


.8e-03 


7 


8. 


194411 


2. 


4e-03 


8 


8. 


194339 


-7 


.le-05 


9 


8. 


198023 


3. 


7e-03 


10 


8. 


200161 


2. 


le-03 


11 


8. 


196456 


-3 


.7e-03 



Table 16: Estimations of the value 5i 



5i tn (uj) = (5 a " n ^'2 ) for the map ([8]) (where S is the 
Feigenbaum constant). The different boxes correspond (from left to right) to E = 10 _1 ,10~ 2 
and 10~ 3 . In all the cases we have taken u = v ^~ 1 . 



• First numerical observation (Section I3.1|) : the sequence a' n (oj) / a' n _ 1 (uj) is not conver- 
gent in re. But, for u) fix, one obtains the same sequence for any family of quasi-periodic 
forced maps like ([9]), with a quasi-periodic forcing of the type g(9,x) = h(x) cos(#). 

• Second numerical observation (Section \3.2\i : the sequence a' n (u)/a' n _ 1 (2uj) is conver- 
gent in n when the quasi-periodic forcing of the type g(9,x) = h(x)cos(9). The limit 
depends on oj and on the particular family considered. 

• Third numerical observation (Section I3.3|) : the two previous observations are not 
true when the quasi-periodic forcing is of the type g(9,x) = gE(0,x) = hi(x)cos(9) + 
Eti2(x) cos (26*) when £ / 0. But the sequence a' n (u>) / a' n _ 1 (2uj) associated to the map ([9]) 
with g = qe is -B-close to the same maps with g = go 



These numerical observations evidence the existence of some structure which govern the asymp- 
totic behavior of the sequence a' n {uj). This structure is indeed the fixed point of a suitable 
renormalization operator acting on the space of functions where the families live. The dynam- 
ics of this operator determine the asymptotic behavior of the sequences a l n (uj)/a' n _ 1 (uj). This 
dynamics depend on the Fourier expansion of the quasi-periodic forcing g(9,x), giving place 
to different behaviors depending on the number of non-trivial Fourier nodes of g(9,x). This is 
described with much more detail in the series of papers |20l l2"Tj |2"2"] . In the first one we give the 
definition of the operator for the case of quasi-periodic maps and we use it to prove the existence 
of reducibility loss bifurcations when the coupling parameter goes to zero. In the second one 
we give a theoretical explanation to each of the numerical observations above in terms of the 
dynamics of the quasi-periodic renormalization operator. Our quasi-periodic extension of the 
renormalization operator is not complete in the sense that several conjectures must be assumed. 
In [22] we include numerical computations which support our conjectures and we show that the 
theoretical results agree completely with the behavior observed numerically. 



References 



[1] D.H. Bailey, Y. Hida, X.S. Li, and Y. Renard. http://crd.lbl.gov/ dhbailey/mpdist/. 



14 



K. Bjerklov. SNA's in the quasi-periodic quadratic family. Comm. Math. Phys., 286(1):137- 
161, 2009. 

E. Castella and A. Jorba. On the vertical families of two-dimensional tori near the triangular 
points of the bicircular problem. Celestial Mech. Dynam. Astronom., 76(l):35-54, 2000. 

W. de Melo and S. van Strien. One-dimensional dynamics, volume 25 of Ergebnisse der 
Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)]. 
Springer- Verlag, Berlin, 1993. 

J. -P. Eckmann and P. Wittwer. A complete proof of the Feigenbaum conjectures. J. Statist. 
Phys., 46(3-4) :455-475, 1987. 

M. J. Feigenbaum. Quantitative universality for a class of nonlinear transformations. J. 
Statist. Phys., 19(l):25-52, 1978. 

M. J. Feigenbaum. The universal metric properties of nonlinear transformations. J. Statist. 
Phys., 21(6):669-706, 1979. 

U. Feudel, S. Kuznetsov, and A. Pikovsky. Strange nonchaotic attractors, volume 56 of 
World Scientific Series on Nonlinear Science. Series A: Monographs and Treatises. World 
Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2006. Dynamics between order and 
chaos in quasiperiodically forced systems. 

J. Figueras and A. Haro. Computer assisted proofs of the existence of fiberwise 
hyperbolic invariant tori in skew products over rotations. Preprint available at 
|http : //www . ma . utexas . edu/mp^arc/c / 10/ 10- 185 . pdf [ 2010. 

C. Grebogi, E. Ott, Pelikan S., and J. A. Yorke. Strange attractors that are not chaotic. 
Phys. D, 13(l-2):261-268, 1984. 

J. F. Heagy and S. M. Hammel. The birth of strange nonchaotic attractors. Phys. D, 
70(1-2):140-153, 1994. 

T. H. Jaeger. Quasiperiodically forced interval maps with negative schwarzian derivative. 
Nonlinearity, 16(4): 1239-1255, 2003. 

A. Jorba, P. Rabassa, and J.C. Tatjer. Period doubling and reducibility in the quasi- 



periodically forced logistic map. Preprint available at http : // www.maia.ub . es/dsg/2011/ 
2011. 

A. Jorba and J. C. Tatjer. A mechanism for the fractalization of invariant curves in quasi- 
periodically forced 1-D maps. Discrete Contin. Dyn. Syst. Ser. B, 10(2-3) :537-567, 2008. 

K. Kaneko. Doubling of torus. Progr. Theoret. Phys., 69(6):1806-1810, 1983. 

O. E. Lanford, III. A computer- assisted proof of the Feigenbaum conjectures. Bull. Amer. 
Math. Soc. (N.S.), 6(3):427-434, 1982. 

M. Lyubich. Feigenbaum-Coullet-Tresser universality and Milnor's hairiness conjecture. 
Ann. of Math. (2), 149(2):319-420, 1999. 

A. Prasad, V. Mehra, and R. Ramaskrishna. Strange nonchaotic attractors in the quasiperi- 
odically forced logistic map. Phys. Rev. E, 57(2): 1576-1584, 1998. 

P. Rabassa. Contribution to the study of perturbations of low dimensional maps. PhD 
thesis, Universitat de Barcelona, 2010. 



15 



[20] P. Rabassa, A. Jorba, and J.C. Tatjer. Towards a renormalization theory for quasi- 
periodically forced one dimensional maps I. Existence of reducibily loss bifurcations. In 
preparation, 2011. 

[21] P. Rabassa, A. Jorba, and J.C. Tatjer. Towards a renormalization theory for quasi- 
periodically forced one dimensional maps II. Asymptotic behavior of reducibility loss bifur- 
cations. In preparation, 2011. 

[22] P. Rabassa, A. Jorba, and J.C. Tatjer. Towards a renormalization theory for quasi- 
periodically forced one dimensional maps III. Numerical support to the conjectures. In 
preparation, 2011. 

[23] D. Sullivan. Bounds, quadratic differentials, and renormalization conjectures. In American 
Mathematical Society centennial publications, Vol. II (Providence, RI, 1988), pages 417-466. 
Amer. Math. Soc, Providence, RI, 1992. 

[24] C. Tresser and P. Coullet. Iterations d'endomorphismes et groupe de renormalisation. C. 
R. Acad. Sci. Paris Ser. A-B, 287(7) :A577-A580, 1978. 

[25] In [20] we address this problem and we prove that this is the case under suitable hypotheses. 

[26] Actually we should say />geometrically asymptotically equivalent, but we will simply say 
asymptotically equivalent. 



16 



